Theoretical simulation of the anisotropic phases of antiferromagnetic thin films 
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We simulate antiferromagnetic thin films. Dipole-dipole and antiferromagnetic exchange inter- 
actions as well as uniaxial and quadrupolar anisotropies are taken into account. Various phases 
unfold as the corresponding parameters, J, D and C, as well as the temperature T and the number 
n of film layers vary. We find (1) how the strength A„ of the anisotropy arising from dipole-dipole 
interactions varies with the number of layers v away from the film's surface, with J and with n; (2) a 
unified phase diagram for all n-layer films and bulk systems; (3) a layer dependent spin reorientation 
(SR) phase in which spins rotate continuously as T, D, C and n vary; (4) that the ratio of the SR to 
the ordering temperature depends (approximately) on n only through (D + A/n)/C, and hardly on 
J; (5) a phase transformation between two different magnetic orderings, in which spin orientations 
may or may not change, for some values of J, by varying n. 
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I. INTRODUCTION 

Magnetic thin films are attracting much interest. 
Some of it derives from applications (in electronics 1 ") 
of ferro/antiferromagnetic layered structures, where bias 
hysteresis arises from interactions at the interfaces of 
filmSi and in nanoparticle layers^ Knowledge of the na- 
ture of the magnetically ordered states, as well as of the 
transitions between them, is important. The spin reori- 
entation (SR) transition is most interesting. Continuous 
SR transitions, in which the direction of the magnetiza- 
tion changes continuously with temperature, were first 
observed in bulk ferrimagnets 4 and in canted spin anti- 
ferromagnets ( AFs) £ Discontinuous SR transitions were 
first discovered in the bulk, in AFs- and in ferromagnets. 7 

Competition of various magnetic anisotropies play de- 
cisive role in SR. Simply put, minimization of the en- 
ergy with respect to direction of the magnetization m 
(or some staggered magnetization m s for AFs) gives the 
physical direction of m at very low temperatures. By 
proper choice of the anisotropy, the energy minimum can 
be controlled, and thus the direction of m. Furthermore, 
in a Ginzburg-Landau like theory, the anisotropy con- 
stants can be made to vary with temperature, and thus 
the direction of m. This approach was first use by Horner 
and Varma&S for continuous SR. Mean field, as well as 
MC calculations, also give continuous SR in the bulkiifii 
Discontinuous SR, on the other hand, occur when one 
local minimum in the free energy, for some spin direc- 
tion, suddenly (as, for instance, the temperature varies) 
becomes the global minimum, at the expense of another 
local minimum, for another spin direction. 

For films, thermally driven SR transitions, whose na- 
ture (continuous or discontinuous) was not clearly es- 
tablished, were first reported for ferromagnets by Pap- 
pas et al»ii Usually?^ but not always^ SR proceeds 
from out of plane to in plane as the temperature in- 
creases. Variation of the number n of layers can also lead 



to SR transitions! 1 ^* 1 ^ 4 Manifestly smooth SR transi- 
tions have been recently observed in ferromagnetic thin 
films There are two main sources for the out of 
plane in-plane anisotropies in films: (1) missing bonds at 
surfaces can give rise to large local magneto-crystalline 
anisotropies then; (2) dipole-dipole interactions induce 
important anisotropies in magnetic films. In ferromag- 
nets, dipolar fields drive spins to lie in plane, rather than 
out of plane, because dipolar field energies (~ m 2 ) that 
obtain when m is out of plane are thus avoided^ 7 " Dipo- 
lar fields lead to stripe like domains in thin films when 
magneto crystalline anisotropies favor spins to be out of 
pkmerl 7 ^!— ^ Growth of such stripes of in plane spins, at 
the expense of out of plane domains (or the other way 
around), as the temperature varies, leads to continuous 
SR in ferromagnetic films m& No continuous SR transition 
is obtained if a homogeneous magnetization as well as 
only a lowest order uniaxial anisotropy are assumedfSiSS 
as has sometimes been done in mean field theoryjSii 2 ^ 
MG 2 ^* 2 ^ and a renormalization group calculation. 26 

The behavior of antiferromagnetic films is qualitatively 
different, mainly because anisotropic effects that arise 
from dipolar fields in AFs are more subtle than in fer- 
romagnets. In AFs, fields decay exponentially beyond 
the system's boundaries, as expected from the following 
simple argument. Consider an AF filling all space where 
z < 0. In the vacuum (i.e., where z > 0), the magnetic 
field h(r) follows from h(r) = V^>(r), where </>(r) is a suit- 
ably defined field. Since 0(r) obeys Laplace's equation for 
z > 0, (f>(r) can be expanded therein, in obvious notation, 
as ^ k akCos(k • r||) exp(— | k | z). This much follows as 
well for ferromagnets. The difference between AFs and 
ferromagnets arises from the fact that whereas Ok — 
for | k \ <\ G | where G gives the periodicity of 4>{x) near 
the surface of an AF, for ferromagnets, | G | scales with 
the inverse ferromagnetic domain size. In addition, the 
previous argument suggests that anisotropic effects from 
dipolar fields may also decrease exponentially, away from 
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surfaces, within AFs. Important qualitative differences 
between anisotropics in ferro- and antifcrromagnets arise 
from this. Unfortunately, relevant experimentsrL22iS£ 
and MC work for one-layer antiferromagnctic films have 
only recently been reportedi22i2i A discontinuous SR has 
been simulated in one layer films with a weak antifer- 
romagnetic exchange^ (in which dipolar interactions are 
dominant) as well as with a strong one.- 1 Because no high 
order (beyond quadratic) site anisotropy was taken into 
account, continuous SRs did not obtain. Finally, there is 
a mean field theory calculation for one-layer Heisenberg 
spin systems which include dipolar interactions as well as 
the lowest order uniaxial anisotropy 3 - which also yields 
a thermally driven discontinuous SR. 

Our aim in this paper is to study (i) how the effective 
surface anisotropy that arises from dipolar interactions 
in magnetically ordered AF films varies with film thick- 
ness and with exchange strength, (ii) how the magnetic 
phases depend on film thickness, as well as on exchange, 
the uniaxial D and quadrupolar C anisotropy constants, 
(iii) how spins on surface layers behave with respect to 
spins on inner layers, and (iv) how the continuous and 
discontinuous SR temperatures depend on various pa- 
rameters. 

The plan of the paper is as follows. The model is spec- 
ified in Sec. [HI Section ITTT1 is about antiferromagnctic or- 
dering in the ground state and the unification that can be 
achieved between film and bulk phase behavior. This uni- 
fication comes about because the anisotropy that arises 
from dipole-dipole interactions is, as surmised in the In- 
troduction, a surface effect. In Sec. IIII Al we define two 
general homogeneous spin configurations. By MC sim- 
ulations, we show that all antifcrromagnetically ordered 
phases, except for the SR phased conform to these con- 
figurations. One (the other one) general configuration 
holds for AF ordered states in which exchange (dipolar) 
interactions dominate. We derive the anisotropy energy 
in each of these two configurations coming from dipolar 
interactions. Monte Carlo results show that the resulting 
effective anisotropy decays exponentially fast with dis- 
tance away from films surfaces. In Sec. IIII Bl the ground 
state continuous SR transition is studied. By MC sim- 
ulations, we study how surface anisotropy arising from 
dipole-dipole interactions drive spin directions as a func- 
tion of layer position. Section IIVI is about thermal ef- 
fects. In Sect. IIV Al we report MC results for transitions 
between various homogeneous magnetic phases. One of 
them is the discontinuous SR transition. In addition, a 
transition between two ordered states, with the same spin 
alignment, is found as the number of film layers changes. 
In Sect. IIVDI we study, by MC simulations, the ther- 
mally driven continuous SR transition. Defining an ef- 
fective uniaxial anisotropy constant D e f f that takes into 
account the dipole-dipole induced anisotropy, we show 
that the ratio of the SR transition temperature to the 
ordering temperature depends on D e ff/C, but depends 
hardly on the exchange constant, as long as it is antifcr- 
romagnetic. 
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FIG. 1: (Color online) (a) Magnetic phase diagram for AFs in 
which exchange dominates. In the SR phase, spins tilt away 
from the z-axis, (f> = ±7r/4, and all spins that are nearest 
neighbors to any one given spin point opposite to it. (b) 
Same as in (a) but for AFs in which exchange is very weak 
and dipolar interactions dominate. A spin configuration in the 
SR phase is depicted in Fig. 2. Full and dashed thick lines 
stand for first- and second-order transitions, respectively. 



II. THE MODEL 

We next specify the model system we study. Let Sj 
be a classical 3-component unit spin at lattice site i of a 
simple cubic (SC) lattice, let 

n = n. r + n d + n A , (1) 

where Hj — ~JJ2{ij) Si ■ Sj, the sum Yltij) ^ s over an 
nearest neighbor bonds, 

w d = E£^f^-> ( 2 ) 

(ij) a/3 
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Yij is the displacement from site i to site j, a is the SC 
lattice parameter, 

Ha = -Dj2(S?f -Cj2i(^) 4 + ( S i)% ( 4 ) 



surface anisotropy that arises from dipole-dipole interac- 
tions in homogeneous states. 



A. Homogeneous states 



and D and C are the uniaxial and quadrupolar anisotropy 
constants, respectively. The nearest neighbor dipolar en- 
ergy Ed is defined through Eqs. (J2J and (j2J. 

The boundary conditions we use are most easily 
grasped in one dimension. Consider first spin sites at 
Xk = ka, for k — — oo . . . , 0, . . . oo. For periodic bound- 
ary conditions (PBC), Sf? = S% +L for all k, and we let 
a spin at the fc-th site interact with all L/2 (L/2 — 1) 
spins immediately to the right (left) of the fc-th site. For 
free boundary conditions (FBC), on the other hand, we 
would let a spin at the A;-th site interact with all spins 
on sites n — 1, . . . , k — 1, k + 1, . . . L. We now return to 
the system of interest here, an n-layer film, by which we 
mean Lx Lxn spins on a fully occupied SC lattice within 
a slab which lies flat on an xy plane. Let the z axis be 
perpendicular to the film layers. We use PBC along the x 
and y directions and FBC along the z direction. Thus, a 
spin on any given site i interacts, through dipolar fields, 
with all other LxLxn-1 spins in the system which are 
in a box, whose top and bottom surfaces coincide with 
the two films surfaces but is otherwise (that is, sidewise) 
centered on the i-th site. 

Our simulations follow the standard Metropolis Monte 
Carlo MC algorithm. 34 In particular, after we choose an 
initial spin configuration, we compute the dipolar field at 
each site. Time evolution takes place as follows. A spin 
is chosen at random and temporarily pointed in a new 
random direction. The move is accepted if either AE < 
0, where AE is the energy change, or with probability 
exp(— AE/ksT), where T is the systems temperature, if 
AE > 0. All dipolar fields are then updated throughout 
the system if the move is accepted, before another spin 
is chosen to repeat the process. By in-plane and out of 
plane we will mean spins lying flat on the xy plane or 
along the z axis, respectively. 



III. EFFECTIVE SURFACE ANISOTROPY 
FROM DIPOLE-DIPOLE INTERACTIONS 

The spin configurations explicitly depicted (not the SR 
phase) in Figs. 1(a) and 1(b) were shown in Ref. [l(| to 
be ground states for large L x L x L systems with PBC. 
Our Monte Carlo calculations show that the same spin 
configurations are also ground states for films, with PBC 
at the film edges and FBC on the perpendicular direction 
to the film. In these states, spins on the two film surfaces 
do not deviate at all from the direction they would point 
to in the bulk. 35 We shall refer to these states as homo- 
geneous. Our Monte Carlo calculations also show (see 
below) that, in the SR phase, spins on surface layers tilt 
away from these directions. We first derive the effective 



Consider first the phases in Fig. QJa), in which, 

S; = (sin 9 cos 4>, sin 9 sin</>, cos 6)rji, (5) 



where 



x(i)+y(i)+z(i) 



(6) 



and x(i), y(i), z(i) is the three dimensional position of the 
i-th site. By proper choice of 9 and <f>, the above equa- 
tions define the three spin configurations shown in Fig. 
1(a). We shall refer to these spin configurations, which 
minimize Ej, as AFj configurations. Spins in these con- 
figurations are clearly collinear. 

Note first that Eqs. © and © imply 



E, 



(3 — —)J, 
n 



(7) 



independently of 9 and <fi, and Ea — —D cos 2 9 — 
C(sin 4 <j> + cos 4 (j>) sin 4 9. To calculate Ed, we substitute 
Eqs. © and © into Eq. © and ©. Note first that 



ViVj = if a 7^ /3 (by x and y reflection symme- 
try). Similarly, TifViVj = Y.,, ' /'"'/<'/..• Therefore 

E d = (2N)- 1 Y^l + (T# - 7 cos 2 9} (8) 



follows, which, by numerical evaluation, gives 

_ 0.67 £d A 2 
E d = cos 9 



where 



n n 



A = 1.984e ri 



(9) 



(10) 



for AFj spin configurations. 

We next calculate E,j + Ed + Ea for the three spin 
configurations shown in Fig. 1(b), given by 

S j = (rf sin 9 cos (j>, rf sin 9 sin (j>, r? cos 9) (11) 

where Tj = [rf , rf , rf] is given by 

Ti = u_\\y{i)+z{i) ; (—iWO+vWi (12) 

We shall refer to the above spin configurations, depicted 
in Fig. lb for some values of 9 and <fi and in Fig. [21 for 
arbitrary 9 and <fi, as AFd configurations. Spins in these 
configurations are in general noncollinear. It is worth 
pointing out that, in L x L x L spin systems in cubic 
lattices, Tld is invariant with respect to both 9 and <f> in 
these AFd configurationsiifli 3 ^ 
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FIG. 2: Spin configuration denned by Eqs. lilt and i I Lit . In 

the SR phase of Fig. 1(b), < 9 < vr/2 and 4> = ±tt/4; for 
other phases of Fig. 1(b), 8 and <f> are as depicted therein. 

Note first that Eqs. {TTJ and imply, 
1 2 J 

Ej = (1 - -)J+ — cos 2 (9, (13) 
n n 

and Ea = —D cos 2 9 — C(sin 4 + cos 4 4>) sin 4 #. We now 
calculate Ed- For systems with complete cubic symme- 
try, that is, with a cubic lattice structure and a cubic 
shape, with the same type of boundary conditions on all 
surfaces, we have showniS that Ed is invariant with re- 
spect to 9 and <fi in Eq. (XH - By the arguments preceding 
Eq. ©, we obtain 

E d = (2N)- 1 ^TpqfTf r/ + m'Tfif - T t frf rf ) cos 2 0] 
y 

(14) 

for n-layer films with PBC at the edges and FBC on 
the top and bottom surfaces. Straightforward numerical 
calculations give 

E d = - (iM -— ) £d + ^^ cos 2 9 (15) 
\ n J n 

for AFd spin configurations. 

Thus, both for AFj and AF d configurations, an effec- 
tive anisotropy, 

D eff = D + A/n (16) 
obtains, where A is given by Eq. 11U|) for AFj, and by 
A = -1.23e d -2J (17) 



for AFd states. 

The anisotropy induced in AFs by dipole-dipole inter- 
actions, given by Eqs. (|10fl . I|16fl . and i|17|) . differ from 
the one for ferromagnets in two respects: (1) it favors 
out of plane spins over in plane ones in AFj states, and 
(2) the corresponding energy for AFs varies as 1/n as 
n increases. The reason for it is given in the Introduc- 
tion, in plane orientation need not be favored in AFs, 
because no significant vacuum dipolar field energy exists 
for them. Wc next discuss the mechanism underlying the 
other effect, the 1/n behavior. 

Equations (|10|l and 117|) suggest that anisotropy effects 
arising from dipole-dipole interactions occur only on sur- 
face layers. We have (numerically) calculated how the 
dipolar field varies with the distance from a film's surface 
for homogeneous spin configurations. Let v = 0, 1, 2, • • • 
number the layers, starting with for one of the two 
outermost surface layers. We find that the deviation of 
dipolar fields from their bulk value decreases exponen- 
tially as v increases. More specifically, 

5h v ~ <5/iexp(— kv), (18) 

where k ~ 4.4 (k ~ 7) for all AFjs and out of plane AFdS 
(for all in plane AFd's) ordered states we have tried. This 
is in agreement with the discussion in Sec. I, since the 
wave vector's magnitude of the field near the surface of 
an AFj or an out of plane AFd (for an in plane AFd) 
on a SC lattice is \f2i:/a (\/5vr/a). Therefore, to the 
accuracy of our numerical results, all anisotropy effects 
arising from dipole-dipole interactions occur only on the 
outer surfaces. Therein, it is given by A/2 for each of 
the two surfaces on n = 2 films, and by A for n = 1 layer 
films. 

In order to characterize antiferromagnetically ordered 
states, we now define 

mj = N- 1 J2S?r ll , (19) 

i 

and 

m a d =N- l Y,STrt. (20) 

i 

In a AFj state, mj = (sin8 cos <f), sin 6* sin </>, cos 9) and 
m d = 0. On the other hand, in a AF d state mj = and 
m d = (sin cos cj>, sin 9 sin <fi, cos 9) . We can define these 
order parameters for the whole system, by summing over 
all sites i, or we can define, say, order parameters for 
surface film layers or interior layers, by summing over 
surface or interior sites. We next discuss the effective 
anisotropy in the SR phase. 

B. The SR phase 

In Ref. the nature of the spin reorientation phase 
[marked as SR in Figs. 1(a) and 1(b)] in L x L x L 
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FIG. 3: (Color online) versus —D e ff/C, for a — x,z on 
two different layers. m^(u) is for the v-th layer (v = for a 
surface layer, and so on). All data follow from MC simulations 
of films of 16 x 16 xn spins, for n = 4, C — — 0.8e,j, and J = 0, 
near T = 0. 



7r/ 2 as D e ff/C sweeps over the 0—1 interval. The order 
parameter, m d and m d , on the surface clearly differs from 
the order parameters on inner layers in Fig. |21 To the 
accuracy of our results, spins on all inner layers do follow 
either Eq. JSJ) or Eq. (|llfl . Thus, phase diagrams for n- 
layer films collapse into a single diagram if D is replaced 
by D + A/n, but only approximately so for the SR phase. 

In order to look further into this effect, we have 
also performed MC simulations of n-layer films with 
anisotropy constants D and C on all sites, except for all 
surface sites, where a variable quantity AD is added to 
D. We calculate | m d {in) | for inner layers, and | m d (s) \ 
for the two surface layers. A plot of | m^iin) — m^s) |, 
for a — x,y, z, vs AD je d is shown in Fig. 0]for T ~ 0, 
D = -e d , C = -2D, and J = 0. | m%{in) - m°t(s) \ 
vanishes for all a at AD = 0.615£<j, as was to be ex- 
pected from Eqs. (|17fl and i|18|) . Further MC simulations 
we have performed for films of various thicknesses yield 
analogous results. 
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AD/8, 

d 

FIG. 4: (Color online) |m^(m) — m^(s)| versus AD for slabs 
of L x L x n dipoles for J = 0, D = —Ed and C = —2D. We 
use FBC in the z direction and PBC in the xy plane. The 
value of D on the two outermost layers of the slab differs by 
AD from the value it takes on the inner layers. 



IV. PHASE TRANSITIONS 

Up to this point we have assumed which of the two, 
AFj or AFd, states the system is in, but we are now able 
to specify which of these two obtain given the value of 
J. Much of this section, the portions having to do with 
phase transitions as the number of film layers change at 
very small temperatures, follow from the following con- 
siderations. Comparison of Eqs. Q, (0, and 11U|) with 
((Toll and O shows that: (1) AFj {AF d ) order ensues 
in out of plane spin configurations when J < — 1.346c; 
(-1.34e d < J < 0); AFj [AF d ] order ensues for in 
plane spin configurations when J < —(1.34 + 0.27/n)ed 
[-(1.34 + 0.27/n)e d < J < 0]. From these conditions on 
| J | one can decide whether | J | is sufficiently small 
for a system to qualify as a dipolar antiferromagnet. 37 
As n — > co, the results obtained in Ref. 01 f° r bulk 
systems follow. Note also that a transformation be- 
tween AFj and AF d ordered states for in-plane config- 
urations can occur as the number n of layers changes if 
-1.61e d <J< -1.34e d . 



systems with PBC was discussed in some detail. In the 
SR phase rrij ^ (m^ ^ 0) for all three components of 
a if J < — 1.34ed (— 1.34^ < J < 0), that is, spins in the 
SR phase are tilted some angle < 8 < ir/2 away from 
the z axis. For homogeneous states, given by either Eq . 
(JSJ or Eq. ip"^ . as for L x L x L systems in Ref. [l0|. 
minimization of the total energy, gives 

S = (±u, ±u, ±v) (21) 

for C < D e ff < 0, where u = y/D e ff/2C and v = 
y/l — 2u 2 . Then, varying D e ff/C over the 0—1 range 
would lead to a n/2 spin rotation. Inspection of Fig. 
[3] shows that the homogeneity assumption is wrong for 
thin films, and, consequently, spins do not quite rotate by 



A. The C=0 transition 

The phase transition at C = for D e ff < is illus- 
trated in Fig. |SJ We know of no previous experimental or 
MC work on this transition in films. It appears to be of 
first order, as predicted by Landau's theory, because no 
symmetry group in any of these phases is a subgroup of 
another one. The transition moves slightly off the C = 
line as T departs from 0, as shown for a one layer film in 
Fig. El giving rise to a reentrant transition. The transi- 
tion at C = 0, however, remains unmoved at T = as the 
number of film layers varies. This is in agreement with 
the statement that dipolar interactions shift the value of 
D but not of C. 
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FIG. 5: (Color online) Order parameter m^, for a = x, y vs C 
for slabs of 16 x 16 x n dipoles on sc lattices at T — 0.05ed/kB , 
where fcs is Boltzmann's constant, for D — —Ed and J = 0. 
A transition between the zy-collinear and xy-canted phases 
is clearly exhibited at C — for different values of n. We 
use FBC in the z direction and PBC in the xy plane. C was 
lowered from C = Ed in AC = — 0.02ed/fcs steps. Each data 
point follows from averages over 10 MC sweeps. 
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FIG. 6: (Color online) Transition temperatures vs —C/D e ff 
from MC simulations of L x L x 1 spin systems, for L = 16 
and L = 32, varying C, with D e ff = — 0.5ed and J — 
fixed. Open and closed symbols stand for L = 16 and L = 32, 
respectively, xy and z stand for in plane and out of plane 
spin configurations. In the simulations, the temperature was 
lowered in steps of 0.02e t j/fcs and O.Oled/ks for L = 32 and 
L = 16, respectively. At each temperature, 10 6 and 2 x 10° 
MC sweeps were taken, for L = 16 and L = 32, respectively. 
In order to check for non-equilibrium effects, the system was 
heated for some values of C / D e ff. No such effects were found. 



B. Discontinuous SR transitions 

Consider first the phase transition between in plane 
and out of plane spin configurations at low temperature 
for C > 0. Assume temporarily the same AFj or AFd 
order in both in plane and out of plane phases, and recall 



that the ground state energy variation with spin tilt angle 
9 is given by Ha if D is replaced by O e // in Eq. 10} . 
Then, 



H A = C 



D, 



If 



C 



D, 



eff' 



(22) 



follows from Eq. replacing D by -D e // and as- 

suming Si = 1 or Si = 1. Clearly, (1) D eff /C < 1 
(D e ff/C > 1) implies an in plane (out of plane) phase in 
the ground state, and (2) the transition is discontinuous 
since Eq. ljL'21) gives an energy barrier between 8 = and 
8 = 7r/2 when D e ff = C. This is as depicted in Figs. 
1(a) and 1(b). It has been observed in experimental and 
numerical work on films The assumption we made, 
that the same AFj or AFd order prevails in both in plane 
and out of plane phases, holds for most values of J, as 
follows from Eqs. Q, Q, ©, ©, D2J), and |Q7J). More 
specifically, the transition occurs between two AFd states 
if —1, 34ed < J < 0. The transition is between two AFj 
states if J < -(1.34 + 0.27/n)e d . Monte Carlo results 
that include this transition are shown in Fig. 0for films 
of n — 1 and n = 2. 

The boundary line between out of plane and in plane 
phases tilts away from D e jf — C as the temperature 
increases, as shown in Fig. \7\ Thus, the possibility 
of thermally driven SR transitions arise, as T varies if 
—0-5£d i; D e ff — C < 0. This is qualitatively as in the 
mean field prediction^ for J = — I0 3 ed, C = 0, and 
n = 1 in Ref. HI. 
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FIG. 7: (Color online) Transition temperatures vs D e ///C 
for n — 1 and n — 2 layer films in which J = — 10, C — 
and C — 1. In the graph, z and xy stand z-collinear and xy- 
collinear, respectively. All data points come from MC simu- 
lations of L x L x n spins for the values of L and n shown 
in the graph. The phase transition boundaries at the top fol- 
low from the location of specific heat peaks obtained while 
lowering the temperature in AT = 0.1 steps. In order to 
make sure equilibrium is realized, the lower phase transition 
boundary, for SR, is obtained from counting, over MC runs 
of several times 10 8 sweeps, the frequency of occurrences of 
the two phases, which must be the same for both phases at 
the boundary line. 
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FIG. 8: (Color online) Everything shown is for J = — 1.38ed, 
C = 0.5£<j, fcsT = O.Oled- This is the lowest temperature 
we obtained with MC simulations in which T was lowered 
in AfcsT = O.lEd steps, down to k B T = 0.01e d - Some 10 5 
MC sweeps were taken at each temperature step, (a) 9 vs D 
for n = 2 (•) and n = 3 (■). The shaded rectangle covers 
values of D/sd where the out of plane AFj state, shown in 
(c) transforms into the the AFd in plane phase, shown in 

(b) if a layer is added to the film, (b) A 9 = 7r/2, AFd, spin 
configuration, on lattice sites of a vertical cut of an n = 3 layer 
film, obtained at fcsT = O.led for C = 0.5 and D = —Q.2ed- 

(c) Same as in (b) but a 9 = 0, AFj, spin configuration that 
obtains for an n = 2 layer film, for the same values of C, D, 
and T, as in (b). 
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FIG. 9: Ground state energies vs 1/n for AFj and AFd states, 
for out-of and in-plane spin alignment, of n-layer slabs, with 
J = — 1.38ed, D = — 0.7ed, and C = 0. • and full line stand 
for in-plane (spins either along the lattice axes or diagonally 
to them) AFd and AFj ordered states, respectively; dashed 
line, and ♦ stand for out of plane AFd and AFj ordered states, 
respectively. Dotted lines are guides to the eye. 



discontinuous SR transitions becomes 



D = C + 2 J + (2.68 - \AA/n)e d , 



(23) 



for C > 0. Note that quantity D e ff is not well defined 
in this narrow J range. We can however use the two dif- 
ferent values D e ff has on both sides of the phase tran- 
sition for comparison of the energies of AFj and AFd 
ordered systems. It can be checked straightfordwardly 
that, again, there is an energy barrier between the 9 = 
and 9 = n/2 phases. This transition, which has not, as 
far as we know, thus far been observed, is illustrated in 
Fig. [SI with MC results for 3 and 2 layer films. The dark 
rectangle shown in Fig. showing the range of values of 
D/sd where n — 2 (n = 3) films order in out of plane AFj 
(in plane AFd) states, follows from our MC simulations. 
It is slightly displaced to the left, by D/ed — 0.12, from 
the prediction that follows from Eq. I|23|l . Irreversibility, 
that keeps spins from reorienting, from in plane to out 
of plane, at low T, as T decreases, is responsible for this 
effect. 



In the small -(1.34 + 0.27/n)e rf < J < -1.34e d range 
the situation is a more interesting. As specified at the top 
of Sect. IIVI the phase transition is then between an out 
of plane AFj ordered state and an in plane AFd ordered 
state. Equation does not apply then, because the 
assumption underlying it, that the same spin ordering 
AFj or AFd prevails on both sides of the phase boundary, 
breaks down. Assuming homogeneity, and making use 
of Eqs. gj), Q), ©, (O, and ||TSJ|, the condition for 



C. Transitions between AFj and AFd ordered states 

Phase transformations that do not involve SR can also 
occur between AFj and AFd ordered states, as the num- 
ber n of layers changes, if —(1.34 + 0.27/n)ed < J < 
-1.34e d . This is illustrated in Fig. Hfor J = -1.38e d , 
D = —0.7ed, and C = 0, where a transition from an 
AFj to an AFd ordered state, both in plane, is shown to 
take place as n decreases from n = 6 to n = 5. (This 
is followed by a spin rotation, from in plane to out of 
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plane, as n decreases from n = 3 to n = 2.) We are 
not aware of any experimental observation of this kind of 
phase transformation. 
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FIG. 10: (Color online) m z d (o and •), m y d (□ and ■), m% (0 

and ♦), and C/fcs (continuous and dashed lines) vs T. Open 
(closed) symbols and continuous (dashed) lines are for MC 
results for systems of 32 x 32 x 1 (of 16 x 16 x 1) spins on sc 
lattices for J — 0, C = — e<j, and D — 0.7(7. At each value 
of T, 4 x 10 5 MC sweeps were made. Lines are guides to the 
eye. 




FIG. 11: Surface and inner layer order parameters | m d (s) | 
and | m2(in) \, respectively, for a — x,y, and z, and specific 
heat, C, vs T. All data points come from MC simulations of of 
8x8x4 spins with J = 0. Symbols (dashed lines) are for order 
parameters (specific heat) of systems with spatially homoge- 
neous anisotropy constants D — —Ed and C = — 2ed and FBC 
on the top and bottom surfaces; continuous lines are for films 
with PBC on all of its boundaries with a uniaxial anisotropy 
constant D — —Sd on inner surfaces, D — —Ed + A/2 on its 
top and bottom surfaces, and C = ~2e d everywhere. 



D. Continuous SR transitions 

A thermally driven SR transition is illustrated in Fig. 
HUI for an n — 1 layer film, ft is rather similar to thermally 
driven transitions in L X L X L systems with PBC^S 
However, as is pointed out in Sec. IIII Bl (and illustrated 
in Fig. EJl for T ~ 0, the SR phase in films with n > 1 
is special. Spin configurations in the SR phase are not 
homogeneous. Whereas spins on inner layers follow Eqs. 
© or l(TT|) . spins on surface layers do not, if D and C are 
homogeneous throughout the system. 

Inhomogeneity effects that arise from the effective sur- 
face anisotropy (induced by dipole-dipole interactions) 
are illustrated in Fig. ^2 for a 4-layer film as a func- 
tion of temperature. Two kinds of data points from MC 
simulations are shown: (1) For films with spatially homo- 
geneous anisotropy constants D and C and FBC on the 
top and bottom surfaces; (2) for films with PBC on all 
of its boundaries 3 ^ with a uniaxial anisotropy constant 
that is A/2 larger on its top and bottom surfaces than 
on the inner layers. Note how the order parameters on a 
surface layer differs from the order parameters on inner 
layers in the former case, and how direct application of an 
anisotropy A/2 on surface layers of films with PBC (no 
anisotropy from dipolar interactions arises then) leads to 
the same effect. This is as expected from the discussion 
in Sec. IIII Bl concerning Fig. 

The continuous SR portion of the phase diagram for 
films is rather like the one for bulk AFsiiS We give the 
MC results we have obtained for n — 1 layer dipolar films 
(in which J = 0) in Fig. [H 



Finally, we compare ordering and SR temperatures. 
Let T z {T xy ) stand for the temperatures below which 
rn z j ^ or m\ ^ (m X j V ^ or m^ y ^ 0). Monte 
Carlo results for the T z /T xy ratio for n = 1 films and 
bulk systems is shown in Fig. El f° r J/^d = and for 
Ed/ J = 0. In the latter case, T z /T xy is, of course, in- 
dependent of J, since J sets the only energy scale then. 
The insensitivity of T z /T xy to dimensionality and to the 
value of J, even when Ed ^ is remarkable. 



T 




FIG. 12: (a) Phases of dipolar (J = 0) antiferromagnetic films 
for C — —Sd- All data points come from MC simulations of 
(o) 16 x 16 x 1 and (x) 32 x 32 x 1 spins. a;y-canted and 
z-collinear stand for the ordered phases that are depicted in 
Fig. 1(b); the SR phase is as pictured in Fig. 2. 
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- D ef/ C 

FIG. 13: (Color online) Temperature ratio T z /T xy vs 
—D e f f/C for the values shown of Ed , J, and C in films (d = 2) 
and Lx Lx L systems with PBC (d — 3). All data points, ex- 
cept the (■) mean field (MF) results, come from MC simual- 
tions. Note that o and A stand for pure dipolar systems 
(J = 0), is for Ed = 0, and □ is for a one layer film in which 
both D and C are constant (D e // varies however, following 
Eqs. 1161 and 1171 . because J varies). Systems are either of 
L x L x L spins (for d=3) or of L x L x 1 spins (for d=2). 
All symbols stand for L = 16 except for o, which stands for 
L = 32. Lines are guides to the eye. 



ial plus quadrupolar anisotropics . We have found how 
the strength A of the effective uniaxial anisotropy that 
arises from dipolar interactions varies with the strength 
of the antiferromagnetic exchange interaction and with 
layer position. We have argued (in Sec. I), and checked 
with MC simulations (in Sec. IIII Ajl . that A decays expo- 
nentially fast with the distance from either of the two out- 
ermost film layers. We have also found (in Sec. IHIj l that, 
except for the SR phase, all antiferromagnetic phases are 
homogeneous, that is, there are no surface states. These 
results, which are peculiar to AFs, imply that all but the 
SR phases of n layer antiferromagnetic films fit into a 
single phase diagram [see Figs. 1(a) and 1(b)] if we let 
D — > D + A/n. The only exception, the SR phase, as- 
sociated with continuous SR (see Sees. IIII Bl and IIV D(l . 
and phase transitions that arise from variation of film 
thickness, in the narrow — 1.61e<2 < J < — 1.3Aed range, 
in which spin orientation may (see Sec. IIV F3|> or may not 
(see Sec. IIV C|) change. Finally, by means of MC simu- 
lations, we have found that the ratio of the continuous 
SR transition to the Neel ordering temperature depends 
very little on J, as illustrated in Fig. ^3 and seems to 
depend on the number n of film layers only through the 
effective uniaxial anisotropy constant D + A/n. 
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